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In the present article we argue that it is possible to find numerical solution of coupled 
magnetic resonance equation for describing wireless energy transmit, as discussed re- 
cently by Karalis (2006) and Kurs et al. (2007). The proposed approach may be found 
useful in order to understand the phenomena of magnetic resonance. Further observa- 
tion is of course recommended in order to refute or verify this proposition. 



1 Introduction 

In recent years there were some new interests in methods 
to transmit energy without wire. While it has been known 
for quite a long-time that this method is possible theoreti- 
cally (since Maxwell and Hertz), until recently only a few 
researchers consider this method seriously. 

For instance, Karalis et al [1] and also Kurs et al. [2] have 
presented these experiments and reported that efficiency of 
this method remains low. A plausible way to solve this prob- 
lem is by better understanding of the mechanism of magnetic 
resonance [3]. 

In the present article we argue that it is possible to find nu- 
merical solution of coupled magnetic resonance equation for 
describing wireless energy transmit, as discussed recently by 
Karalis (2006) and Kurs et al. (2007). The proposed approach 
may be found useful in order to understand the phenomena of 
magnetic resonance. 

Nonetheless, further observation is of course recommend- 
ed in order to refute or verify this proposition. 

2 Numerical solution of coupled-magnetic resonance 

equation 

Recently, Kurs et al. [2] argue that it is possible to repre- 
sent the physical system behind wireless energy transmit us- 
ing coupled-mode theory, as follows: 

am(t) = {iUm - r m) a m (t ) + 

+ ^ ^ IK'nmQ'n (^) ^m(^) • 

n^m 



These equations can be expressed as linear 1st order ODE 



as follows: 


/'(*) = -iaf(t) +m g(t) 


(4) 


and 


g‘{t) = — iPg(t) 


(5) 


where 


a = {ui i - iTj) 


(6) 


and 


p = (lj 2 - ir 2 ) 


(7) 



Numerical solution of these coupled-ODE equations can 
be found using Maxima [4] as follows. First we find test when 
parameters (6) and (7) are set up to be 1. The solution is: 



(%i5) ’diff(f(x),x)+%i*f=%i*b*g(x); 

(%o5) ’diff(f(x),x,l)+%i*f=%i*b*g(x) 

(%i6) ’diff(g(x),x)+%i*g=%i*b*f(x); 

(%o6) ’diff(g(x),x,l)+%i*g=%i*b*f(x) 

(%i7) desolve([%o5,%o6],[f(x),g(x)]); 

The solutions for /( x) and g(x) are: 

,, , [ig{0)b-if{x)]sm{bx) 

/M = j 

[g(x)~ f(0)b]cos(bx) g(x) 
b + —’ (8) 

, , [if(0)b - ig(x)]sm(bx) 

s(P = j 

_ [/O) -ff(0)fr] cos (bx) + f(x) 



The simplified version of equation (1) for the system of 
two resonant objects is given by Karalis et al. [1, p. 2]: 



Translated back to our equations (2) and (3), the solutions 
for a = P = 1 are given by: 



and 



da i 
dt 
da 2 
dt 



i (cj i — zTi) di + i(ca 2 , 
i (. tu 2 - iV 2 ) a 2 + inai . 



( 2 ) 

(3) 



ai(f) = 



[ia 2 ( 0 )k — ia i] sin(«;f) 



[a 2 -a 1 (0)Ac]cos(/cf) | q 2 

K K 
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f(x) = e 



— ( ic—ia)t/2 



[2z/(0) c + 2ig(0)b - /(0)(zc - za)] sin ^ Vc 2 2 ac+Ab^+a^ ^ 



■J c 2 — 2 ac + 4 6 2 + i 



+ 



+ 



5(a;) = e 



—(ic—ia)t / 2 



/(0) COS ( Vc2 — 2ac+4ft2 +a2 ^ 

■J c? — 2ac + 4 fr 2 + a 2 
[2z/(0) c + 2ig(0) a — g(0)(ic — za)] sin f Vc 2 ^2ac+46 2 +a 2 ^ 



■J c 2 — 2 ac + 4 b 2 + a 2 

5 ( 0 ) cos ( Vc‘i-2ac + 4b'i+ al 
+ ■ 



i Jc 2 — 2 ac + 4fo 2 4- a 2 



a x (t) = e 



-{ip-i a )t /2 ( [2zai(0)/3 -f 2za 2 (0) k — (z/3 — za) a x ] sin(|f) _ ai(0) cos (ft) 
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a 2 (f) = e 



_ -{i/3-ia 



)z /2 /" [ 2 * 02 ( 0 ) p + 2ia 1 (0)K - (ip - za)a 2 ] sin(|f) a 2 (0) cos(| f) 









(13) 



(14) 

(15) 

(16) 



and 

[zoi(0)k — za 2 ] sin(/ct) 

/■C 

[a x - a 2 (0)/c] cos(/cf) + a x 

K K 

Now we will find numerical solution of equations (4) and 
(5) when a + P + 1. Using Maxima [4], we find: 

(%il2) ’diff(f(t),t)+%i*a*f(t)=%i*b*g(t); 

(%ol2) ’diff(f(t),t,l)+%i*a*f(t)=%i*b*g(t) 

(%il3) ’diff(g(t),t)+%i*c*g(t)=%i*b*f(t); 

(%ol3) ’diff(g(t),t,l)+%i*c*g(t)=%i*b*f(t) 

(%il4) desolve([%ol2,%ol3],[f(t),g(t)]); 

and the solution is found to be quite complicated: these are 
formulae (13) and (14). 

Translated back these results into our equations (2) and 
(3), the solutions are given by (15) and (16), where we can 
define a new “ratio”: 

i = \]P 2 - 2 aP + 4/t 2 + a 2 . (12) 

It is perhaps quite interesting to remark here that there is 
no “distance” effect in these equations. 

Nonetheless, further observation is of course recommend- 
ed in order to refute or verify this proposition. 
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